Adiabatic Condition and Quantum Geometric Potential 
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In this paper, we present a (/(l)-invariant expansion theory of the adiabatic process. As its application, we 
propose and discuss new sufficient adiabatic approximation conditions. In the new conditions, we find a new 
invariant quantity referred as quantum geometric potential (QGP) contained in all time-dependent processes. 
Furthermore, we also give detailed discussion and analysis on the properties and effects of QGP. 
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Since the establishment of the quantum adiabatic theorem 
in S [110] in 1923, many fundamental results have been ob- 
tained, such as Landau-Zener transition Q], the Gell-Mann- 
Low theorem 0], Berry phase Qj and holonomy [[s]]. Also 
the adiabatic processes find their applications in the quantum 
control and quantum computation I.9..10. .1 1..12j] . Recently the 
common-used quantitative adiabatic condition Ill5l[l6l[l7[| has 
been found not able to guarantee the validity of the adiabatic 
approximation [13> JJO- Consequently various new conditions 
are conjectured and a series of confusions and debates arise. 
For example, it was argued lisll that the traditional adiabatic 
condition did not have any problem at all and that the invali- 
dation of the condition did not mean the invalidation of adia- 
batic theorem fl9tl . Some new conditions proposed in 1I20II21I1 
but too rigorous to be used conveniently. Although [22] also 
adopted the adiabatic perturbation expansion but did not give 
out proper condition because the basis in 1I22I1 can not show 
certain geometric properties in the adiabatic process. I23I1 
pointed out the limitation of traditional condition but also did 
not give out a proper condition. To solve the problem of in- 
sufficiency of traditional adiabatic condition in fl3i [Till and 
clarify the subsequent confusions, in this paper, we present 
two new sufficient conditions in which the properties and ef- 
fects of a new invariant quantity are detailedly discussed. 

Let us consider a quantum system governed by a time de- 
pendent Hamiltonian H{t) and the initial state of the sys- 
tem is an eigenstate |m, 0) of H(Q) with eigenvalue £,,,(0), 
where m denotes the initial value of dimensionless quantum 
number set. By introducing a dimensionless time parame- 
ter T - E,„ (0) t/H and a dimensionless Hamiltonian Ii{t) = 
H{t)/E,„ (0), the time dependent Schrodinger equation reads 



r = /2(t)|0„,(t)), |0„,(0)) = \m, 0) . 

OT 



(1) 



The exact solution \<^m (t)) to Eq.© is referred to as the sys- 
tem's dynamic evolution orbit in the Hilbert space. 



Furthermore, by considering r as a fixed parameter, we can 
always solve the following quasi-stationary equation 



h (t) \(f„ (t)) = e„ (r) \(p„ (r)) . 



(2) 



And the eigenstate |^„(t)) with the corresponding initial state 
|«, 0) is referred to as the adiabatic orbit of the system. 

For convenience, we denote y„„, = i{ipn{T) \ i^,„(t)) and the 
dot here and below expresses the derivative with respect to 
time. Apparently, an adiabatic orbit multiplied by an arbitrary 
time-dependent phase factor still describes the same adiabatic 
orbit. It is not difficult to see that the following adiabatic orbit 



)) = exp{-/]| 



{e,„{A) - y„„n{A)]dA |(^„,(t)> (3) 



is invariant under the following t/(l) transformation 

\ip,n(T)) ^ e'A'f^' |^„,(t)> (/„(0) = 0). (4) 

Here f,„(Q) = is because of given initial state. We call 
this adiabatic orbit with special choice of the time-dependent 
phase factor as the U{l)-invariant adiabatic orbit. 

It is clear that, although the initial conditions \m, 0) are the 
same, the dynamic evolution orbit |<1),„(t)) do not always coin- 
cide with the adiabatic orbit |i/j„,(t)), or they are not even close 
to each other Obviously they coincide if and only if 



7mn = (V« 9t m). 



(5) 
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Generally speaking, the dynamic evolution orbit \^,„ (t)) 
starting from the initial state |m, 0) will change among some 
adiabatic orbits which will cause transitions between different 
orbits. Our task is to find the proper condition under which the 
dynamic orbit is sufficiently close to the adiabatic orbit when 
the Eq.® is not satisfied. Since the Hamiltonian /2(t) is Her- 
mitian, all the U{1)- invariant adiabatic orbits in Eq.lO at a 
given time constitute a complete orthonormal basis of the sys- 
tem. In this basis, the dynamic evolution orbit of the system 
reads 

|1>™ (t)> = Yj ^"(^^ l*^"'"" (^^) ' 1*^"' = I'"' 0> ■ 
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The expansion in Eq.® is referred to as the U(l)-invariant 
adiabatic expansion with the time-dependent coefficients. 
Therefore, the set of coefficients equations reads 



(7) 



where the diagonal elements of the matrix M(t) are zero and 
the non-diagonal elements of M(r) read 

M(t)„,„ = / ((Df (t) I cbf (r)) = |7,„„(T)| (8) 

where 

Omn{r) = I dA {em(A) - e„(A) + y„„ - rmm)+argy„„,(T). (9) 
Jo 

From Eq.(|7]), it is not hard to get 

C(t)("+i> = A(t)C(t) 



where 



A(t) = 



Z U/,+ l)=«+l p= 

/;=1 



(10) 



(11) 



and (« + 1) denotes the (n -i- 1)''' derivation of (?(t) with respect 
to time. Then we have following theorem 

Theorem For an N-level quantum system and an arbitrary 
real s and a time period T, if the following conditions hold 



maxA(0),j and maxA(r),j < B < oo 



(12) 



iTnuXT")! can be expanded as 2 fl* .e"^™"'" 
with max a*,,, - D < oo 



'ini:m,k 



max 

k I 



2pN(N-\)BD 



= s' < I with s + s' < I. (14) 



then the probability of finding dynamical orbit in the adiabatic 
orbit \^"J'(t)^ is greater than (1 - with 6 = e/(l - e')- 
Proof: From Eq.(|7]i and the conditions above, we have 



k„,(r)-i| = 



2 drLir) i ai„e'-'"-- 



P 



2 2 '^mn 2 ( 1)'' !'■ i ^'+ 

iii=mk=l g=0 
p oo 



c„(0)"" 



(15) 



Namely, 



i-|c,„(r)|<|i-c„,(r)|< 



1 -s' 



(16) 



Therefore, the probability of finding dynamical orbit in the 
adiabatic orbit |<l>^"''(r)^ is 



(17) 



Thus we prove the theorem. 

Although Eq.( 12- 14) in the theorem are sufficient, however, 
it is somewhat too complicated. It is not hard to find for 
any A^-level Hamiltonian with both time-independent terms 
of \jmn\ and 9i„„ satisfying Rydberg-Ritz Combination Prin- 
ciple(RRCP) 0ni + dim - Qnm, for an arbitrary real < 5 <sc 1, 
when the following condition holds 



max 

Vm and k,ni^m 



\ykn 



VIZ., 



\ykm\ 



Vm and Kni^m \e„{T) - e,„(T) + A„„(t)| - 1 



(18) 



(19) 



where 

d 

Km (t) = Tmm (t) - Jnn (t) + arg j„„, (t) (V« ^ m). (20) 

dT 

then the probability of finding dynamical orbit in the adiabatic 
orbit |<l);|f'(T)^ is greater than (1 - 5)^. 

Proof: Denote c^,(t) - e"^'"'^c„i{T) with w„, - a)„ - Omn- 
Then, c^jC") satisfy equations 



i— C'(t) = nc'(T) 

OT 



(21) 



where C'(t) - {c[{t),c^{t), . . . ,c'j^(t))^ , IT is a self-adjoint 
matrix, II^ji = (±>k and 11^:; - lyt/l. Denote eigenvalues of n as 
T]m, we have 12511 



\Tln 



■ C0,„\ < 



(22) 



If unitary matrix U diagonalizes 11, then UTIU^ 
diag{rix,ri2, - ■ ■ ,riN], that is Uik^kj - rjiUij, thus Uij 
2 —Uik. When condition ([Till holds, then 



Md^'ii + i (23) 



?7/ - VA^- 1 V - 1 

so \Uij\ <[6+ V2o(5-)) / and \Uii\ > 1 - 5/2 - o(5-). 

C'{t) can be solved exactly as 11 is time-independent 



C'(t) = e'"^C'(0). 



(24) 



Applying initial condition cj,(0) = 5jt„,, the exact solution of 
<,(''") is 2 It^mil^e""'- Thus, |c,„(t)| > 1 -5, then the probability 

A- 

of finding dynamical orbit in the adiabatic orbit |<I>5jf'(T")) is 

Prni^) > (1 - Thus the proof is completed. 

The premises of Eq.(fT9]l on Hamiltonian are non-trivial. 
For any general 2D system h{T) - e+(T)|+, t)(+, t| + 
e_(T)|-, t)(-, t|, after applying a transformation t ^ t' = 
^"'(t) with ^(t) oc ^ I (+, A\ jj |-, A) \dA, then we can forcibly 
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get e'_ — e'^ + A'^_ = constant which makes the final time- 
dependent 2D system satisfying the Umitations of Eq.(fT9]l. 

In condition Eq.(fT9]l, there appears a new interesting quan- 
tity A,„„ referred to as quantum geometric potential (QGP) 
for following three reasons. First, QGP is also L'(l)-invariant 
under the transformation Eq.(|4li. Second, the integral of QGP 
over a closed smooth curve is the difference of Berry phases of 
different adiabatic orbits. And the last reason is that the value 
of QGP depends only on the path and measure of adiabatic or- 
bit or, in other words, A,„„(r)/ \ymn\ (Vn m) is invariant un- 
der any transformation t — > r' = /(t). Furthermore, It can be 
proved that in 2D systems A„,„(t)/2 |y,„„| is just the geodesic 
curvature of spherical curve corresponding to the adiabatic or- 
bit on the surface of Bloch sphere or 2D real Ray space. 

Proof: Generally, we can write the Hamiltonian of a 2D 
system as /i(t) - A{t) + B(t)«(t) ■ cr, where n(T) - 
(sin 6 (t) cos ip (r) , sin 9 (t) sin if (t) , cos 9 (r)). Choosing ap- 
propriate phases, the Hamiltonian's instantaneous eigenstates 
or adiabatic orbits read 



+, r) = cos ^ |0> + e'^^^' sin 2^ 
-, t) = sin ^ |0) - e'f^^^ cos ^ 



0(t) 



(25) 



It's quite clear that polarization vectors of the above two adi- 
abatic orbits point to «(t) and -n(T) at time t, respectively. 
Considering the adiabatic orbit \+,t), the QGP of this orbit 
can be easily calculated as 



) sin + 29 (p cos 9 + 4> sin" 9 cos 9 - 4>9 sin S 



+ ((/isine) 



(26) 



As a comparison, we will calculate the geodesic curvature of 
the spherical curve r (t) - n{T). 

1^ dr\ d~r 



i sin -H 29~(h cos -H (i^ sin^ 9 cos ( 



Hm9 



^6*2 + ^0 sin 0)^1 



■,(27) 



where curve element 



ds - \dr^ - sji) + {^(p sin fl) dr - 2\y„„\dT. 



Then we get 



(28) 



In the following part, two models will be presented to show 
Eq.(fT9]) is a good sufficient adiabatic condition and the effect 
of QGP is significant. Firstly, we shall study a spin-half parti- 
cle in a magnetic field. The Hamiltonian of the system is 

h(T) = Tja-, + ^ [o-., cos (2Kr]T) + (Ty sin {2Kr]T)\ , (29) 

where rf - hu{)IE±, ^ - fiCL)/E+ and Ej- - 
^jrp- + ^~ are all constants. Obviously, Eq.(fT9l) is a sufficient 



adiabatic condition for this kind of Hamiltonian. Properly 
choosing phases, two adiabatic orbits can be written as 



|^4r)> = cos(f)|0) + e2'-^''^sin(f)|l) 



|^_(T)) = sin(f)|0)-e2'^''^cos(f)|l> 



(30) 



where cos0 = 77/ -y/^^TjT^ Consider the adiabatic orbit 
|^+(t)), we have QGP, A+_ = 2Kricos9. It is easy to obtain 
the expression of the new adiabatic condition of Eq.(fT9ll 



rp- + ^ - Krj cos 9: 



» \KT]sin9\ . 



(31) 



Suppose the initial state of the system is |+,0), we have the 
fidelity between the dynamic evolution orbit and the adiabatic 
orbits at time t 

F(t) = ^cos2(At) + sin^(AT) ^ i^-^M^^^ (32) 



where A = -^(1 - K)^ rf- + is also a constant parameter 

If we choose ?/ » ^ and K ^ \, then the traditional con- 
dition [16] is satisfied but the new condition Eq.(fT9l) is not. 

Meanwhile, the fidelity F {t) ^ Vl ~ cos^ 9 sin" (Ar) -+> 1 
when T is not too small. Thus, even though the traditional con- 
dition is satisfied and we might regard the system as slowly 
changing one, the quantum adiabatic approximation may be 
unfaithful description of the system because of the effect of 
the QGP 

While if we choose ?? » ^ with K ^ \ and » //, in this 
case, the QGP is much larger than the difference of the instan- 
taneous energy eigenvalues, and the new condition Eq.(fT9l) 
is satisfied while the traditio nal one is not. Now we have 
F(t) * Vl -sin^e sin" (At) * 1. Therefore, the QGP can 
help to guarantee the validity of the adiabatic approximation 
despite the difference of energy eigenvalues is too small to 
satisfy the traditional condition. 

Next, notifying that if QGP has same sign as the corre- 
sponding difference of energy eigenvalue, it will positively 
guarantee the system evolution to be adiabatic. Moreover, 
if |A„„,| / K« I m>| » 1, with Eq.(fT9l), then the evolution 
may be adiabatic whether the adiabatic orbit moves slowly 
or fast. Thus we shall present an interesting Hamiltonian 
for illustrating QGP may be helpful to construct robust sys- 
tem. Consider a 2D system governed by Hamiltonian h{T) - 
T]cr, + e-''^'^'^{rjQa^ + T/ie'''^'^' V^e^^'-'^'Oe"'"^'^ with 770/77 » 
1 and 770/771 » 1. The density matrixes of adiabatic orbits 
read 



Pf (r) 4 ± ^ 



■ e-"''^''^(77oo-,. + riie''^^''^^a^e-'''^"''^Y'"^^^ 

2N{t) ' 
(33) 

where A^(t) - ^t]q + (rj + rj] cos 2t]2t)^ + rj^ sin" 2772r. The 
density matrixes of the evolution orbits starting from the cor- 
responding initial states of adiabatic orbits reads 



l\ ^TT, Ji , '7oO"a- + (?7i +?7)o-z '|,,T, ^ 

P±(t) = —U(t) 1 + U (t), 

H±\ ) 2 '\ N{Q) ' 



(34) 
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where U(t) = e-'';'^cTg,>,2<r,Tg-,((;,o+^2)<r,+7,nr,)T ^j^^j fj^g jj^j^^j 
states are p'^'{Q). The probabiHties of staying in the corre- 
sponding adiabatic orbits are 

P±iT) = 2^7(^77(7) (''0+ '71+'?^ cos 2772T+ 277/71 COS^T/lT 

277"* . T _ 477(770 + 772)77- 2 • 2 - 

+ sin 77T — cos 772T sin 77T 

77"^ 77"^ 

77n^ - \ 1 

H — ^ sin 277T sin 27727" + -. (35) 



Here 77 = ^J']l + im + m)^ and 77 = V77770 + 77772+772771. 
Since 770/77 » 1 and 770/771 » 1, then probabiHties will ob- 
tain a lower bound f independent on the magnitude of 772: 
- 1 - (77 + 771)^ /N(0)-, which approach to 1. It's not hard 
to verify when 772 >> 77, A+_ has the same sign as - E+, 
and |A+_| / 1(+ | -)| ^ 770/771 » 1. When 772 is large, the ve- 
locity of the adiabatic orbit has the same order of magnitude 
of 772, at this time, the adiabatic orbit fast oscillates around the 
exact dynamic evolution orbit, but the evolution of the system 
still keeps adiabatic. Fig. 1 shows evolution orbit and adiabatic 
orbit for 770/77 = 770/771 = 20 and 770/772 - 0.2. 




For a short summary, it is worthwhile to point that, by the 
theorem or new adiabatic condition Eq.(fT9l), the problems 
showed in [13,14] has not existed because the relation be- 
tween systems a and b constructed in [13,14] does not guar- 
antee them. Apparently, different from those conditions in 
ll2ll 122! I23I1 . our conditions are presented in a more natural 
way full of geometric interpretation. One more hint we may 
get here is that we should more carefully deal with the phase 
appearing in the time-dependent evolution. It is just improp- 
erly handling the phase of Eq.® in the work of predecessors 
ifTsi [l6ll that led to their improper traditional condition and 
later contradiction in iflsifTill . The condition ( fT9] l also implies 
a modification of the difference of energy eigenvalues is nec- 
essary. Description of the time-dependent evolution might be 
more precise and more appropriate via replacing e„, (T)-e„ (t) 
by e,„ (t) - e„ (t) + A„,„. And a related experiment ll24ll for 
verifying the effect of QGP has been finished. The experi- 
ment also found the characteristic frequency of a kind of time- 
dependent systems should be corrected via QGP. The experi- 
ment also illustrated the QGP should reflect some properties 
of time-dependent systems, and is not just a convenient mathe- 
matical technique. As it is shown in our paper, QGP may play 
an important role in some kinds of time-dependent procedure, 
but what role it may play in general time-dependent system 
is not clear now. We guess non-trivial QGP will more or less 
affect the evolution procedure of time-dependent system. 

In conclusion, according to the concepts of [/(l)-invariant 
adiabatic orbit and t/(l) invariant expansion stated in this pa- 
per, we present a theorem and a new sufficient adiabatic con- 
dition, from which we get an interesting quantity QGP with 
its effects and geometric properties detailedly discussed. At 
the end we present two models to show the significant effect 
of QGP on the evolution. 



FIG. 1: evolution orbit(red line) and adiabatic orbit(blue line). 

This kind of models allows the parameters of the system 
have a certain variant range, as long as the adiabatic condi- 
tion Eq.(fT9]l holds. Thus we may conclude that the QGP may 
help setting up robust systems which may tolerate faults of the 
system itself. Another interesting hint from this model is that 
the adiabatic orbit may be very complicated comparing with 
evolution orbit, which is counterintuitive from the traditional 
opinion. 
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